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We report new limits on the ty° — > jX decay of the neutral pion into a photon and a new gauge 
boson X followed by the decay X — > e + e~ . If this process exist, one would expect a flux of high 
energy X's produced from 7r°'s generated by the proton beam in a neutrino target. The X's would 
then penetrate the downstream shielding and be observed in a neutrino detector via their decays. 
Using bounds from the NOMAD and PS191 neutrino experiments at CERN that searched for an 
excess of e + e~ pairs from heavy neutrino decays, stringent limits on the branching ratio as small as 
Br(w° -> jX) < 10" 16 are obtained. These limits are several orders of magnitude smaller than the 
previous experimental and cosmological bounds. The obtained results are used to constrain models, 
where the X interacts with quarks and leptons, or it is a new vector boson mixing with photons, that 
transmits interaction between our world and hidden sectors consisting of SU(3)c x SU(2)l X U(1)y 
singlet fields. 

PACS numbers: 14.80.-j, 12.60.-i, 13.20.Cz, 13.35.Hb 



Many extensions of the Standard Model (SM) such as 
GUTs super-symmetric |2|, super-string models [3[ 
and models including a new long-range interaction, i.e. 
the fifth force predict an extra U (1) factor and there- 
fore the existence of a new gauge boson X correspond- 
ing to this new group. The predictions for the mass of 
the X boson are not very firm and it could be light 
enough (Mx <C Mz) for searches at low energies. If 
the mass Mx is of the order of the pion mass, an effec- 
tive search could be conducted for this new vector boson 
in the radiative decays of neutral pscudoscalar mesons 
P — > jX, where P — ir°, r\, or rj' , because the decay rate 
of P — >• 7 + any new particles with spin or -| proves 
to be negligibly small Q. Therefore, a positive result in 
the direct search for these decay modes could be inter- 
preted unambiguously as the discovery of a new light spin 
1 particle, in contrast with other experiments searching 
for light weakly interacting particles in rare K, tt or fi 
decays 0-0] - 

The decay ir° — > 7 A can also occur in several interest- 
ing extensions of the SM that suggest the existence of, so- 
called, "hidden" sectors consisting of SU (3)o x SU (2) L x 
U(1)y singlet fields. These are sectors of weakly interact- 
ing massive particles (WIMPs), that do not interact with 
the ordinary matter directly and couple to it by gravity 
and possibly by other very weak forces. If the mass scale 
of a hidden sector is too high, it will be experimentally 
unobservable. However, there is a class of models where 
Dark Matter WIMPs have interactions with the SM me- 
diated by Uh(l) gauge group and the corresponding hid- 
den gauge boson could be light, see e.g. |8rfl0|. The 
fact that the X coupling strength may be in the experi- 
mentally accessible and theoretically interesting regions, 
makes further searches for X's very attractive at the high 
intensity frontier [TTI |. 

From the analysis of data from earlier experiments, 
constraints on the branching ratio for the decay of P — >• 
7 + X range from 10 -7 to 10 -3 depending on whether 



X interacts with both quarks and leptons or only with 
quarks [12j. Direct searches for a signal from ir° — > 7 A" 
decay have been performed in a few experiments with 
two different methods: i) searching for a peak in inclu- 
sive photon spectra from two-body tt° — > 7 + nothing 
decays, where "nothing" means that X is not detected 
because it either has a long life time or decays into vv 
pairs [l3l - [l7| . and ii) searching for a peak in the invari- 
ant mass spectrum of e + e~ pairs from ir° decays, which 
corresponds to the decay X — > e + e~ (lH [l7|. 

The best experimental limit on the branchning ratio of 
the decay tt° ->• jX, Br(TT Q -> jX) < (3.3 - 1.9) x 1CT 5 
(90%C.£.) for X masses ranging from to 120 MeV, 
was obtained by the NOMAD experiment at CERN [lj]. 
Using 450 GeV proton collisions with the SPS neutrino 
target as a source of high energy X's from 7r° — > jX 
decays, they searched for a signal from the Primakoff 
conversion X — > n° in their detector [l9| . The best bound 
for the branching ratio Br(ir° — > r )X)Br(X — > e + e~) < 
2 x 10~ 4 - 4 x 10~ 6 for the X mass range 25 < M x < 
120 MeV was obtained by the SINDRUM collaboration 
at PSI [l6|. In this letter we show that more stringent 
limits on the decay ir° — > -fX, followed by the decay 
X — > e + e~ can be obtained from the results of sensitive 
searches for an excess of single isolated e + e~ pairs from 
decays of heavy neutrinos in the sub-GeV mass range by 
the PS191 [H [2l| and NOMAD [H-H experiments at 
CERN. 

Consider first the NOMAD experiment on search for 
the decay i>h ve^e~ of a heavy neutrino i>h into a 
lighter neutrino and e + e~ pair performed at the CERN 
West Area Neutrino Facility (WANF) 0. The WANF 
provides an essentially pure beam for neutrino exper- 
iments. It consists of a Be target irradiated by 450 GeV 
protons from the CERN SPS. The secondary hadrons are 
focused by the horn and the reflector, and protons that 
have not interacted in the target, secondary hadrons and 
muons that do not decay are absorbed by a 400 m thick 
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shielding made of iron and earth. The NOMAD detec- 
tor described in Ref. [25[ is located at 835 m from the 
neutrino target. It consists of a number of sub-detectors 
most of which are located inside a 0.4 T dipole magnet 
with a volume of 3.5 x 3.5 x 7.5 m 3 including an active 
target of drift chambers (DC) with a mass of 2.7 tons. 
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FIG. 1: Schematic illustration of a proton beam dump exper- 
iment on search for tt° — > jX — » e + e~ decay chain: neutral 
pions generated by the proton beam in the neutrino target 
(T) produce a flux of high energy X's which penetrate the 
downstream shielding and decay into e + e~ pair in a neutrino 
detector. 

If the decay ir° — > 7X exists, one expects a flux of high 
energy X bosons from the WANF target, since 7r°'s are 
abundantly produced in the forward direction by high 
energy protons either in the target or in the beam dump 
following the decay tunnel. If A is a relatively long- 
lived particle, this flux would penetrate the downstream 
shielding without significant attenuation and would be 
observed in the NOMAD detector via the X — > e + e~ 
decay into a high energy e + e _ pair, as schematically il- 
lustrated in Fig. [TJ The occurrence of X — >• e + e~ de- 
cays would appear as an excess of isolated e + e _ pairs in 
NOMAD above those expected from standard neutrino 
interactions. The experimental signature of these events 
is clean and they can be selected with small background 
due to the excellent NOMAD capability for precise mea- 
surements of e + e~ pairs, see example [3, H^- As the 
final states of the decays Vh — > ve + e~ and X — > e + e~ 
are identical, the results of the searches for the former 
can be used to constrain the later for the same e + e _ 
invariant mass regions. 

The calculated flux and energy spectrum of ir° pro- 
duced in the Be target were used to predict the flux of 
X's with momenta pointing to the NOMAD fiducial area 
as a function of the X mass, see Fig. [2] The search for 
the Vh ve + e~~ decay described in |22J corresponds to 
the total number of 4.1 x 10 19 protons on target (pot). 
The strategy of the analysis was to identify — > ^e + e~ 
candidates by reconstructing in the DC isolated low in- 
variant mass e + e~ pairs that are accompanied by no 
other activity in the detector. The measured rate of 
e + e~ pairs was then compared to that expected from 
known sources. By selecting only two tracks identified as 
an e + e~ pair with energy > 4 GeV and invariant mass 
M e + e - < 95 MeV, the inverse total momentum pointing 
back to the target, and forming a vertex within the DC 
fiducial volume, an excess of 2.1 events compatible with 
background expectations was observed. For a given flux 
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FIG. 2: Combined energy spectrum of X bosons with mass 
Mx = 10 MeV from the SPS neutrino target and from the 
beam dump region at the NOMAD detector calculated for 
Br(Tv° jX) = 1. 



d$(Mx , Ex , Np t) jdEx of A's the expected number of 
X — > e + e~ decays occuring within the fiducial length L 
of the NOMAD detector located at a distance L' from 
the neutrino target is given by 
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where Ex , Px , and tx are the X energy, momentum 
and the lifetime at rest, respectively, e is the e + e~ pair 
reconstruction efficiency, N pot is the number of pot. The 
acceptance A of the detector was calculated tracing X's 
produced in the Be-target or beam dump to the detector 
taking the relevant momentum and angular distributions 
into account. As an example for a mass Mx = 10 MeV, 
A = 8.5% and e ~ 20%. The X flux calculated as a 
function of the X energy is shown in Fig. [2] 

The final 90% C.L. exclusion region in the Br(n° — > 
jX)Br(X — > e + e~) vs tx plane shown in Fig. [3] to- 
gether with the PS191 result (see below) is calculated 



by using the relation N^+^L > N x ^ e +e-i where N' 
(= 2.1 events) is the 90% C.L. upper limit for the ex- 
pected number of signal events [22| ■ Our result is sensi- 
tive to a branching ratio Br(n° — > jX) > 10~ 15 , which 
is more than nine orders of magnitude smaller than the 
previous limit from the experiment (l6j . Over most of 
the tx region, the X lifetime is sufficiently long, that 
LMx/pxtx < L'Mx/pxtx < 1. 

To obtain limits for the region Mx > 95 MeV, we con- 
sider next the results of the PS191 neutrino experiments 
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FIG. 3: The 90% C.L. upper limits on the branching ratio 
Br(-K° -» jX)Br(X -s> e+e") versus r x from the NOMAD 
(solid) and PS191 (dashed, dotted) experiments. The num- 
bers near the curves indicate the corresponding values of Mx ■ 



at CERN that searched for decays Vh — > ve + e of heavy 
neutrinos produced in 2-body decays tt, K — >■ e, 11 + Vh in 
the v h mass range from ~ 10 to ~ 350 MeV [13 This 
experiment, specifically designed to search for neutrino 
decays in a low-energy neutrino beam, was performed 
by using 19.2 GeV protons from the CERN Proton Syn- 
chrotron with the total number of 0.86 x 10 19 pot. The 
PS191 detector, located at the distance of 128 m from the 
target, consist of a 12 m long decay volume, eight cham- 
bers located inside the volume to detect charged tracks 
and followed by a calorimeter. The events searched for 
in the experiment were requested to consist of two tracks 
originating from a common vertex in the decay volume 
and giving rise to at least one shower in the calorimeter. 
No single e + e~ events were observed and limits were es- 
tablished on the v e{i — Vh mixing strength as a function 
of the Vh mass. 

Using Eq.([T|) for the PS191 case, similar to above con- 
siderations, we calculte the 90 % C.L. branching ratio 
limit curve shown in Fig. [3] together with the result 
from NOMAD. For the mass region 95 < M x < 135 
MeV/c 2 the best limits from PS191 are in the region 
Br(7T° 1 X)Br(X -> e+e") < (2 - 4) x 10~ 15 . In 
this estimate the average X momentum is < px >— 1 
GeV and the decay region length is I = 7 m. The results 
obtained were also cross checked with the limits set by 
PS191 on the v e-fl — Vh mixing strength [2(| and found to 
be in agreement. For the mass region below 95 MeV, NO- 
MAD provides better bounds than PS191. The n° -)• jX 
decay rate constrained by NOMAD is small enough to 
produce a detectable excess of e + e~ events in the PS191 
experiment. In addition, for masses Mx < 10 MeV the 



reconstruction efficiency drops, because the opening an- 
gle of the e + e _ pair is ~ Mx/Ex < 10~ 2 rad, which 
is small enough to be resolved in the PS191 detector. 
Hence, some events would be misidentified as a single 
track and, would be rejected. 

The obtained results can be used to impose constraints 
on the magnitude of the coupling constant ax = g\ 1^ 
for the interaction of X bosons with both quarks and 
leptons, which can be written in the form: 

L x = 9x{QbxB 1 + Q eX Ll + Q^xLl + Q T xL l T )X l (2) 

where B % = E 9 =«,d, s ,.. 97*9, K = e^e + v eLl l v eLl 
see e.g. [HQ- Assuming charges Qbx — Qex — 1, we 
found 



ax < 3.4 x 10 



-13 



V Ml J 



M x [MeV] 



-3/2 



(3) 



which is valid for Mx < M w , where M^ is the mass of 
tt , and t x > 10~ 10 M x [MeV] s. This limit is more re- 
strictive than those obtained in @, 0| j a- n d than bounds 
reported by NOMAD [l8j]. Less stringent limits (by a 
factor ~ a) could be obtained for the cases where the X 
interacts only with leptons, or when it is a leptophobic 
boson which interacts only with quarks and decays dom- 
inantly into e + e _ pair through the quark loop if its mass 
M x > 2m e [10. 

We can also constrain recent models, where new 
vector-like bosons mediate interaction between our world 
and hidden sectors due to the kinetic mixing term 



L int = -\ X F^. 



(4) 



where, \ is the mixing strength, and F^ v , V^ v are the 
ordinary photon and hidden vector field strengths, re- 
spectively, see e.g. |10| . As follows from ((4]), any source 
of 7's could produce kinematically possible massive states 
V according to the appropriate mixings. If the mass of V 
is below the mass of 7r°, it can be produced in the decay 
tt° — > jX, with the subsequent decay of X into e + e~ 



pair 



The corresponding branching ratio is given by: 



Br(ir° jV) = 2x 2 (l 



MI ) 
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For V masses below the 7r°-meson mass, My < M„, the 
dominant V boson decay is e + e~ with a rate which, for 
small mixing, is given by: 
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Using Eqs.(5,6) we can determine the 90% C.L. exclusion 
regions in the (My; y) plane from the results of the NO- 
MAD [22j and PS191 [2(| experiments, which are shown 
in Fig. @] together with the area exclude d by the electron 
beam dump experiment E137 at SLAC [28|, see Ref-Hi] 



4 




10 

M v , MeV 

FIG. 4: 90% C.L. exclusion regions in the (My, \) plane 
obtained from the results of the NOMAD PS191 [H, 
and positronium [3(j experiments. The areas excluded by the 
electron beam dump experiment E137 [28|] and by the KLOE 
experiment are also shown for comparison. 



and discussion therein, positronium experiment search- 
ing for the decay e + e _ — > 7 V [3(3] and recent results from 
KLOE 31]. The shape of the exclusion contour from the 
PS 191 experiment corresponding to the X mass range 
My > 95 MeV is defined mainly by the phase space 
factor in ([5]). One can see, that the NOMAD exclusion 
region falls basically within that of the E137 experiment. 
However, the fact that the NOMAD events are generated 
by the proton beam, while the E137 events are originated 
from the electron beam dump is important if, e.g. the X 
is a leptophobic boson which interacts dominantly with 
quarks. The attenuation of the X- flux due to X inter- 
actions with matter was found to be negligible, e.g. for 
couplings of (j3|) the X boson mean free path in iron is 3> 



100 km, as compared with the iron and earth shielding 
total length of 0.4 km used for the NOMAD beam. 

In summary, using sensitive limits from the PS191 and 
NOMAD experiments on heavy neutrino decays Uh ^ 
ve + e~ in the sub-GeV mass range, new stringent limits 
on the branching fraction Br(ir° — > jX) < 10~ 8 - 10~ 15 , 
for X masses 1 < M x < 135 MeV and lifetime 10~ n < 
tx < 1 s, are obtained. Our best result is sensitive to a 
branching ratio Br(ir° —> jX) > 10 -15 , which is about 
nine orders of magnitude smaller than the previous limit 
from the SINDRUM experiment jl6j]. It is also a factor 
~ 100 more stringent than the bound from cosmologi- 
cal considerations (321 • The obtained results are used to 
set new limit on the X coupling strength to lepton and 
quarks, and also to constrain models, in which mixing 
between photons and a new vector-like bosons mediates 
interaction between our world and hidden sectors con- 
sisting of SU(3) C x SU(2) L x U(1) Y singlet fields. For 
the X mass range < M w the obtained limits on the mix- 
ing x are more stronger than those recently derived by 
the experiment KLOE [3l|, and comparable with those 
obtained in Ref. [29l l34j . Stringent constraints on mixing 
strength 77 obtained from SM1987A observations, have 
been recently reported in 33j. We demonstrate a sig- 
nificant potential of experiments similar to NOMAD for 
the future more sensitive searches for new physics at 
the high intensity frontier [111 ]. For example, stringent 
limits for the extended sub-GeV X boson mass range 
Mx > M n could be obtained from the searches for 77, 77' 
mesons decays 77, 77' — > jX with the subsequent decays 
X —> e + e~, 7r + 7r~, // ± e T , /i ± 7r =F , ... We note that 

limits on the decays 77, 77' — > jX could also be obtained 
from PS191 or other existing data provided cross-sections 
for 77, 77' productions in the forward direction in p — Be 
collisions are known. These limits are expected to be 
more restrictive than those recently reported by the Al 
collaboration from the search for the decay X — > e + e~ 
in a fixed target experiment [35j |. 

I am indebted to N.V. Krasnikov for useful discussions, 
and M.M. Kirsanov for help in calculations. 
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